3.415 \(\int \cos ^4(c+d x) (a+a \sec (c+d x)) (A+B \sec (c+d x)+C \sec ^2(c+d x)) \, dx\)

Optimal. Leaf size=102 \[ \frac{a (A+B+C) \sin (c+d x)}{d}+\frac{a (3 A+4 (B+C)) \sin (c+d x) \cos (c+d x)}{8 d}-\frac{a (A+B) \sin ^3(c+d x)}{3 d}+\frac{1}{8} a x (3 A+4 (B+C))+\frac{a A \sin (c+d x) \cos ^3(c+d x)}{4 d} \]

[Out]

(a*(3*A + 4*(B + C))*x)/8 + (a*(A + B + C)*Sin[c + d*x])/d + (a*(3*A + 4*(B + C))*Cos[c + d*x]*Sin[c + d*x])/(
8*d) + (a*A*Cos[c + d*x]^3*Sin[c + d*x])/(4*d) - (a*(A + B)*Sin[c + d*x]^3)/(3*d)

________________________________________________________________________________________

Rubi [A]  time = 0.210859, antiderivative size = 102, normalized size of antiderivative = 1., number of steps used = 7, number of rules used = 6, integrand size = 39, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.154, Rules used = {4074, 4047, 2635, 8, 4044, 3013} \[ \frac{a (A+B+C) \sin (c+d x)}{d}+\frac{a (3 A+4 (B+C)) \sin (c+d x) \cos (c+d x)}{8 d}-\frac{a (A+B) \sin ^3(c+d x)}{3 d}+\frac{1}{8} a x (3 A+4 (B+C))+\frac{a A \sin (c+d x) \cos ^3(c+d x)}{4 d} \]

Antiderivative was successfully verified.

[In]

Int[Cos[c + d*x]^4*(a + a*Sec[c + d*x])*(A + B*Sec[c + d*x] + C*Sec[c + d*x]^2),x]

[Out]

(a*(3*A + 4*(B + C))*x)/8 + (a*(A + B + C)*Sin[c + d*x])/d + (a*(3*A + 4*(B + C))*Cos[c + d*x]*Sin[c + d*x])/(
8*d) + (a*A*Cos[c + d*x]^3*Sin[c + d*x])/(4*d) - (a*(A + B)*Sin[c + d*x]^3)/(3*d)

Rule 4074

Int[((A_.) + csc[(e_.) + (f_.)*(x_)]*(B_.) + csc[(e_.) + (f_.)*(x_)]^2*(C_.))*(csc[(e_.) + (f_.)*(x_)]*(d_.))^
(n_)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)), x_Symbol] :> Simp[(A*a*Cot[e + f*x]*(d*Csc[e + f*x])^n)/(f*n), x]
+ Dist[1/(d*n), Int[(d*Csc[e + f*x])^(n + 1)*Simp[n*(B*a + A*b) + (n*(a*C + B*b) + A*a*(n + 1))*Csc[e + f*x]
+ b*C*n*Csc[e + f*x]^2, x], x], x] /; FreeQ[{a, b, d, e, f, A, B, C}, x] && LtQ[n, -1]

Rule 4047

Int[(csc[(e_.) + (f_.)*(x_)]*(b_.))^(m_.)*((A_.) + csc[(e_.) + (f_.)*(x_)]*(B_.) + csc[(e_.) + (f_.)*(x_)]^2*(
C_.)), x_Symbol] :> Dist[B/b, Int[(b*Csc[e + f*x])^(m + 1), x], x] + Int[(b*Csc[e + f*x])^m*(A + C*Csc[e + f*x
]^2), x] /; FreeQ[{b, e, f, A, B, C, m}, x]

Rule 2635

Int[((b_.)*sin[(c_.) + (d_.)*(x_)])^(n_), x_Symbol] :> -Simp[(b*Cos[c + d*x]*(b*Sin[c + d*x])^(n - 1))/(d*n),
x] + Dist[(b^2*(n - 1))/n, Int[(b*Sin[c + d*x])^(n - 2), x], x] /; FreeQ[{b, c, d}, x] && GtQ[n, 1] && Integer
Q[2*n]

Rule 8

Int[a_, x_Symbol] :> Simp[a*x, x] /; FreeQ[a, x]

Rule 4044

Int[csc[(e_.) + (f_.)*(x_)]^(m_.)*(csc[(e_.) + (f_.)*(x_)]^2*(C_.) + (A_)), x_Symbol] :> Int[(C + A*Sin[e + f*
x]^2)/Sin[e + f*x]^(m + 2), x] /; FreeQ[{e, f, A, C}, x] && NeQ[C*m + A*(m + 1), 0] && ILtQ[(m + 1)/2, 0]

Rule 3013

Int[sin[(e_.) + (f_.)*(x_)]^(m_.)*((A_) + (C_.)*sin[(e_.) + (f_.)*(x_)]^2), x_Symbol] :> -Dist[f^(-1), Subst[I
nt[(1 - x^2)^((m - 1)/2)*(A + C - C*x^2), x], x, Cos[e + f*x]], x] /; FreeQ[{e, f, A, C}, x] && IGtQ[(m + 1)/2
, 0]

Rubi steps

\begin{align*} \int \cos ^4(c+d x) (a+a \sec (c+d x)) \left (A+B \sec (c+d x)+C \sec ^2(c+d x)\right ) \, dx &=\frac{a A \cos ^3(c+d x) \sin (c+d x)}{4 d}-\frac{1}{4} \int \cos ^3(c+d x) \left (-4 a (A+B)-a (3 A+4 (B+C)) \sec (c+d x)-4 a C \sec ^2(c+d x)\right ) \, dx\\ &=\frac{a A \cos ^3(c+d x) \sin (c+d x)}{4 d}-\frac{1}{4} \int \cos ^3(c+d x) \left (-4 a (A+B)-4 a C \sec ^2(c+d x)\right ) \, dx+\frac{1}{4} (a (3 A+4 (B+C))) \int \cos ^2(c+d x) \, dx\\ &=\frac{a (3 A+4 (B+C)) \cos (c+d x) \sin (c+d x)}{8 d}+\frac{a A \cos ^3(c+d x) \sin (c+d x)}{4 d}-\frac{1}{4} \int \cos (c+d x) \left (-4 a C-4 a (A+B) \cos ^2(c+d x)\right ) \, dx+\frac{1}{8} (a (3 A+4 (B+C))) \int 1 \, dx\\ &=\frac{1}{8} a (3 A+4 (B+C)) x+\frac{a (3 A+4 (B+C)) \cos (c+d x) \sin (c+d x)}{8 d}+\frac{a A \cos ^3(c+d x) \sin (c+d x)}{4 d}+\frac{\operatorname{Subst}\left (\int \left (-4 a (A+B)-4 a C+4 a (A+B) x^2\right ) \, dx,x,-\sin (c+d x)\right )}{4 d}\\ &=\frac{1}{8} a (3 A+4 (B+C)) x+\frac{a (A+B+C) \sin (c+d x)}{d}+\frac{a (3 A+4 (B+C)) \cos (c+d x) \sin (c+d x)}{8 d}+\frac{a A \cos ^3(c+d x) \sin (c+d x)}{4 d}-\frac{a (A+B) \sin ^3(c+d x)}{3 d}\\ \end{align*}

Mathematica [A]  time = 0.403058, size = 97, normalized size = 0.95 \[ \frac{a (24 (3 A+3 B+4 C) \sin (c+d x)+24 (A+B+C) \sin (2 (c+d x))+8 A \sin (3 (c+d x))+3 A \sin (4 (c+d x))+24 A c+36 A d x+8 B \sin (3 (c+d x))+48 B c+48 B d x+48 C d x)}{96 d} \]

Antiderivative was successfully verified.

[In]

Integrate[Cos[c + d*x]^4*(a + a*Sec[c + d*x])*(A + B*Sec[c + d*x] + C*Sec[c + d*x]^2),x]

[Out]

(a*(24*A*c + 48*B*c + 36*A*d*x + 48*B*d*x + 48*C*d*x + 24*(3*A + 3*B + 4*C)*Sin[c + d*x] + 24*(A + B + C)*Sin[
2*(c + d*x)] + 8*A*Sin[3*(c + d*x)] + 8*B*Sin[3*(c + d*x)] + 3*A*Sin[4*(c + d*x)]))/(96*d)

________________________________________________________________________________________

Maple [A]  time = 0.109, size = 141, normalized size = 1.4 \begin{align*}{\frac{1}{d} \left ( Aa \left ({\frac{\sin \left ( dx+c \right ) }{4} \left ( \left ( \cos \left ( dx+c \right ) \right ) ^{3}+{\frac{3\,\cos \left ( dx+c \right ) }{2}} \right ) }+{\frac{3\,dx}{8}}+{\frac{3\,c}{8}} \right ) +{\frac{Aa \left ( 2+ \left ( \cos \left ( dx+c \right ) \right ) ^{2} \right ) \sin \left ( dx+c \right ) }{3}}+{\frac{Ba \left ( 2+ \left ( \cos \left ( dx+c \right ) \right ) ^{2} \right ) \sin \left ( dx+c \right ) }{3}}+Ba \left ({\frac{\cos \left ( dx+c \right ) \sin \left ( dx+c \right ) }{2}}+{\frac{dx}{2}}+{\frac{c}{2}} \right ) +aC \left ({\frac{\cos \left ( dx+c \right ) \sin \left ( dx+c \right ) }{2}}+{\frac{dx}{2}}+{\frac{c}{2}} \right ) +aC\sin \left ( dx+c \right ) \right ) } \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(cos(d*x+c)^4*(a+a*sec(d*x+c))*(A+B*sec(d*x+c)+C*sec(d*x+c)^2),x)

[Out]

1/d*(A*a*(1/4*(cos(d*x+c)^3+3/2*cos(d*x+c))*sin(d*x+c)+3/8*d*x+3/8*c)+1/3*A*a*(2+cos(d*x+c)^2)*sin(d*x+c)+1/3*
B*a*(2+cos(d*x+c)^2)*sin(d*x+c)+B*a*(1/2*cos(d*x+c)*sin(d*x+c)+1/2*d*x+1/2*c)+a*C*(1/2*cos(d*x+c)*sin(d*x+c)+1
/2*d*x+1/2*c)+a*C*sin(d*x+c))

________________________________________________________________________________________

Maxima [A]  time = 0.937492, size = 178, normalized size = 1.75 \begin{align*} -\frac{32 \,{\left (\sin \left (d x + c\right )^{3} - 3 \, \sin \left (d x + c\right )\right )} A a - 3 \,{\left (12 \, d x + 12 \, c + \sin \left (4 \, d x + 4 \, c\right ) + 8 \, \sin \left (2 \, d x + 2 \, c\right )\right )} A a + 32 \,{\left (\sin \left (d x + c\right )^{3} - 3 \, \sin \left (d x + c\right )\right )} B a - 24 \,{\left (2 \, d x + 2 \, c + \sin \left (2 \, d x + 2 \, c\right )\right )} B a - 24 \,{\left (2 \, d x + 2 \, c + \sin \left (2 \, d x + 2 \, c\right )\right )} C a - 96 \, C a \sin \left (d x + c\right )}{96 \, d} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(cos(d*x+c)^4*(a+a*sec(d*x+c))*(A+B*sec(d*x+c)+C*sec(d*x+c)^2),x, algorithm="maxima")

[Out]

-1/96*(32*(sin(d*x + c)^3 - 3*sin(d*x + c))*A*a - 3*(12*d*x + 12*c + sin(4*d*x + 4*c) + 8*sin(2*d*x + 2*c))*A*
a + 32*(sin(d*x + c)^3 - 3*sin(d*x + c))*B*a - 24*(2*d*x + 2*c + sin(2*d*x + 2*c))*B*a - 24*(2*d*x + 2*c + sin
(2*d*x + 2*c))*C*a - 96*C*a*sin(d*x + c))/d

________________________________________________________________________________________

Fricas [A]  time = 0.504472, size = 221, normalized size = 2.17 \begin{align*} \frac{3 \,{\left (3 \, A + 4 \, B + 4 \, C\right )} a d x +{\left (6 \, A a \cos \left (d x + c\right )^{3} + 8 \,{\left (A + B\right )} a \cos \left (d x + c\right )^{2} + 3 \,{\left (3 \, A + 4 \, B + 4 \, C\right )} a \cos \left (d x + c\right ) + 8 \,{\left (2 \, A + 2 \, B + 3 \, C\right )} a\right )} \sin \left (d x + c\right )}{24 \, d} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(cos(d*x+c)^4*(a+a*sec(d*x+c))*(A+B*sec(d*x+c)+C*sec(d*x+c)^2),x, algorithm="fricas")

[Out]

1/24*(3*(3*A + 4*B + 4*C)*a*d*x + (6*A*a*cos(d*x + c)^3 + 8*(A + B)*a*cos(d*x + c)^2 + 3*(3*A + 4*B + 4*C)*a*c
os(d*x + c) + 8*(2*A + 2*B + 3*C)*a)*sin(d*x + c))/d

________________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(cos(d*x+c)**4*(a+a*sec(d*x+c))*(A+B*sec(d*x+c)+C*sec(d*x+c)**2),x)

[Out]

Timed out

________________________________________________________________________________________

Giac [B]  time = 1.24325, size = 294, normalized size = 2.88 \begin{align*} \frac{3 \,{\left (3 \, A a + 4 \, B a + 4 \, C a\right )}{\left (d x + c\right )} + \frac{2 \,{\left (9 \, A a \tan \left (\frac{1}{2} \, d x + \frac{1}{2} \, c\right )^{7} + 12 \, B a \tan \left (\frac{1}{2} \, d x + \frac{1}{2} \, c\right )^{7} + 12 \, C a \tan \left (\frac{1}{2} \, d x + \frac{1}{2} \, c\right )^{7} + 49 \, A a \tan \left (\frac{1}{2} \, d x + \frac{1}{2} \, c\right )^{5} + 28 \, B a \tan \left (\frac{1}{2} \, d x + \frac{1}{2} \, c\right )^{5} + 60 \, C a \tan \left (\frac{1}{2} \, d x + \frac{1}{2} \, c\right )^{5} + 31 \, A a \tan \left (\frac{1}{2} \, d x + \frac{1}{2} \, c\right )^{3} + 52 \, B a \tan \left (\frac{1}{2} \, d x + \frac{1}{2} \, c\right )^{3} + 84 \, C a \tan \left (\frac{1}{2} \, d x + \frac{1}{2} \, c\right )^{3} + 39 \, A a \tan \left (\frac{1}{2} \, d x + \frac{1}{2} \, c\right ) + 36 \, B a \tan \left (\frac{1}{2} \, d x + \frac{1}{2} \, c\right ) + 36 \, C a \tan \left (\frac{1}{2} \, d x + \frac{1}{2} \, c\right )\right )}}{{\left (\tan \left (\frac{1}{2} \, d x + \frac{1}{2} \, c\right )^{2} + 1\right )}^{4}}}{24 \, d} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(cos(d*x+c)^4*(a+a*sec(d*x+c))*(A+B*sec(d*x+c)+C*sec(d*x+c)^2),x, algorithm="giac")

[Out]

1/24*(3*(3*A*a + 4*B*a + 4*C*a)*(d*x + c) + 2*(9*A*a*tan(1/2*d*x + 1/2*c)^7 + 12*B*a*tan(1/2*d*x + 1/2*c)^7 +
12*C*a*tan(1/2*d*x + 1/2*c)^7 + 49*A*a*tan(1/2*d*x + 1/2*c)^5 + 28*B*a*tan(1/2*d*x + 1/2*c)^5 + 60*C*a*tan(1/2
*d*x + 1/2*c)^5 + 31*A*a*tan(1/2*d*x + 1/2*c)^3 + 52*B*a*tan(1/2*d*x + 1/2*c)^3 + 84*C*a*tan(1/2*d*x + 1/2*c)^
3 + 39*A*a*tan(1/2*d*x + 1/2*c) + 36*B*a*tan(1/2*d*x + 1/2*c) + 36*C*a*tan(1/2*d*x + 1/2*c))/(tan(1/2*d*x + 1/
2*c)^2 + 1)^4)/d